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The higher-dimensional vacuum Einstein equation admits translationally non-uniform 
black string solutions. It has been argued that infinitesimally non-uniform black strings 
should be unstable in 13 or fewer dimensions and otherwise stable. We construct nu- 
merically non-uniform black string solutions in 11, 12, 13, 14 and 15 dimensions. Their 
stability is investigated using local Penrose inequalities. Weakly non-uniform solutions 
behave as expected. However, in 12 and 13 dimensions, strongly non-uniform solutions 
appear to be stable and can have greater horizon area than a uniform string of the same 
mass. In 14 and 15 dimensions all non-uniform black strings appear to be stable. 

1 Introduction 

Consider vacuum General Relativity in D dimensions compactified on a Kaluza-Klein 
(KK) circle. This theory admits a black string solution given by the product of a D — 1 
dimensional Schwarzschild solution of radius r + with a circle of circumference L = 2irR. This 
black string has horizon topology S D ~ 3 x S 1 and is invariant under translations around the 
S 1 . Gregory and Laflamme showed that the string is unstable against linearised gravitational 
perturbations if r + < r* for some critical radius r* proportional to R For example, if 
D = 5 then r* = 0.8762.R. 

Some time ago, it was conjectured that this theory also admits non-uniform black strings 
(NUBSs), which break the translational invariance around S 1 , and form a 1-parameter family 
which bifurcates from the uniform black string family at the marginally stable solution with 
critical radius r* [2j. This conjecture has been confirmed by perturbative construction of 
such solutions, valid for small non-uniformity [31 H], and by numerical solution of the Einstein 
equation [SJ El El- 

A convenient parameterization of these solutions is with the non-uniformity parameter [3] 
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where R max and R m i n are the maximum and minimum values of the radius of the S D ~ 3 of the 
horizon. It is expected that non-uniform black strings exist for all A > 0, and that at A = oo, 
there is a merger between the family of non-uniform black strings, and a 1-parameter family 
describing black holes localized on the KK circle, with horizon topology S D ~ 2 [8]. Numerical 
evidence in support of this expectation has been obtained for D = 5, 6 [9j [TDl U\- 

The perturbative construction reveals an interesting change in the behaviour of non- 
uniform black strings at the critical dimension D = 13 [I]. For D < 13, infinitesimally 
non-uniform black strings have a smaller horizon area than a uniform black string of the same 
mass. This implies that such a solution cannot be the endpoint of the Gregory-Laflamme in- 
stability of the uniform string. However, for D > 13, an infinitesimally non-uniform solution 
has a greater horizon area than a uniform string of the same mass, and is a natural candidate 
for the endpoint of the GL instability. 

The aim of this paper is to investigate the stability of non-uniform black strings with 
finite A. The perturbative results just discussed suggest that weakly non-uniform (i.e. small 
A) black strings should be unstable for D < 13 and stable for D > 13. The stability of 
strongly non-uniform black strings has not been investigated but it has been suggested that 
such strings might be stable for D < 13 [TT] . 

Before we can investigate stability, we must determine the relevant non-uniform black 
string solutions numerically. We are especially interested in strongly non-uniform strings 
with D ~ 13. Such solutions have not been constructed previously. Ref. [6] gave results for 
6 < D < 11. We will construct solutions with D > 11. For D — 11 our solutions have much 
greater A than those of Ref. [U] . 

Our solutions exhibit some novel features. For D < 13, all non-uniform strings constructed 
previously have the property that, as A increases, the mass increases. However, for D = 12, 13 
we find that the mass increases to a maximum at A = Ai and then decreases for A > AiQ The 
first law implies that the horizon area exhibits the same qualitative behaviour as the mass, 
i.e., it is maximized at A = Ai- Furthermore, there exists A2 > Ai such that a NUBS with 
A > A2 has greater horizon area than a uniform string of the same mass. Hence such NUBSs 
are thermodynamically preferred in the microcanonical ensemble. 

Figure [l] gives a sketch of these results]^] For D = 12, 13, the occurrence of a maximum 
mass/area along the NUBS branch results in a cusp when horizon area is plotted against mass. 
We indicate also the expected behaviour of the (not yet constructed) branch of localized black 
hole solutions. 

The behaviour of the NUBS branch that we have discovered for D = 12, 13 is similar to the 
behaviour of localized black hole solutions for D = 5,6. Refs. [TO], [7j found that as one moves 
along the branch of such solutions, starting from a small black hole, the mass increases to a 
maximum and then decreases before the merger with the NUBS branch. The resulting phase 
diagram is shown in Figure [2j with the cusp appearing on the localized black hole branch. 

It would be very interesting if a turning point (i.e. a maximum of mass/area) exists along 

1 Ref. |12j found an analogous feature in D = 5, 6. However, subsequent work [6] could not reproduce this 
feature and neither can we. We believe that the result of Ref. [12] arises from numerical error. 

2 We have done the plots in Figures [l] [2] and [3] with the real data (for the NUBS and the uniform black 
strings) but the different phases are so close to each other that they cannot be easily distinguished by eye. 
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Figure 1: Schematic plot of horizon area against mass for black holes and black strings with 
D = 12, 13. The dashed curve is the uniform black string branch, the solid curve the non- 
uniform string branch and the dotted curve the localized black hole branch. We find that a 
maximum mass/area solution occurs on the non-uniform black string branch, resulting in a 
cusp on the black curve. The localized black hole curve is conjectural but it seems likely that 
this will merge with the NUBS curve. 

the NUBS branch for D < 11. However, our results for D = 11 suggest that this is not the 
case. Our solutions extend all the way up to A ~ 2.4 and we have not seen any evidence 
of a turning point. In fact, when we plot the physical quantities against A, we find that for 
A > 1 they appear to be nearly constant, which is the expected behaviour as one approaches 
a merger point. Identifying the value of -Do such that solutions with D < D have the same 
qualitative behaviour as for D = 5, 6 (i.e. with a turning point along the localized black 
hole branch instead of the NUBS branch) will require constructing the localized black hole 
solutions, but our results show D < 12 and suggest that D = 11. 

Turning to the case of D > 13, we have obtained NUBS solutions for D = 14, 15. This is 
the first time such solutions have been constructed non-perturbatively. We do not discover 
any surprises with these solutions: as A increases, the mass decreases with no sign of any 
turning point. This suggests that the phase diagram should take the form of Figure [3j as has 
been suggested previously [6]. 

We will investigate stability using the method of local Penrose inequalities of Ref. [13], 
which we will review below. This method reveals that a necessary condition for stability of 
a static black hole is that it should be a local maximum of horizon area at fixed mass in the 
space of initial data for Einstein's equation obeying the relevant asymptotic boundary condi- 
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Figure 2: Same plot as figure [T] but for D = 5,6, based on results of Refs [TQl (ZJ- The 
maximum mass solution (the cusp) occurs on the localized black hole branch. Our results 
suggest that this will be the qualitative behaviour for all D < 11. 

tions. Therefore an instability can be demonstrated by constructing initial data describing a 
perturbed black hole which has greater horizon area than an unperturbed black of the same 
mass. 

Using local Penrose inequalities, we find that D < 11 NUBSs are unstable for all values 
of A for which we have constructed solutions. This suggests that all D < 11 NUBSs are 
classically unstable. 

For D = 12, 13, we find that all solutions with A < Ai are unstable. For A > Ai, we 
have constructed multi-parameter families of initial data and found them all to respect the 
local Penrose inequality, which strongly suggests that such solutions are stable]^] Hence for 
D = 12, 13, there is a range of masses M for which there exist stable non-uniform black string 
solutions. In fact, there can exist two stable black string solutions: one uniform and one 
non- uniform with the non-uniform string entropically preferred for A > A2. 

For our D = 14, 15 NUBSs, we have investigated multi-parameter families of initial data 
and found them all to be consistent with the local Penrose inequality, for all values of A for 
which we have constructed solutions. This strongly suggests that these solutions are stable. 

The organisation of the rest of the paper is as follows. In §2 we will describe our numerical 
construction of NUBSs in arbitrary number of dimensions. Our construction is based on Ref. 

3 An instability for A close to, but less than, A x is suggested by the more heuristic "turning point method", 
which has been used to predict neutron star instabilities [2]. For black holes, the local Penrose inequality 
explains the success of this method. 
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Figure 3: Same plot as figure [T] but for D > 13. The behaviour of the localized black hole 
branch is conjectural. 



[7J, which we briefly review in £2.1 . In £2.2 we present our results for the different values of 
D. In £j3]we apply the local Penrose inequality method of [T5] to NUBSs. The results on the 
stability of these solutions are presented in §3.2 



2 Numerical construction of NUBSs 
2.1 Method 

In this section we outline our numerical construction of NUBSs in an arbitrary number of 
space-time dimensions D. We will follow the numerical method of [7J based on the harmonic 
formulation of the Einstein equations for finding static black hole space times (see Chapter 
10 in [UJ for a review). 

Following |7J we write the most general ansatz for the metric with the desired isometries 
(U(1) T x SO(D - 2) in our case)f] 

2 S A 2 A B 

ds 2 = 4 rl A y 2 e T dr 2 + d{l 2 D _ 3) + e A dx 2 + T ° (dy + y f(y) F dx) 2 , (2) 

where = 1 — £ 2 and A = jjjhgf an d ?"o are constants. Here the functions X = 

{T, S, A, B, F} depend on the coordinates (x, y) and they are our unknowns. Setting T = S = 



4 Here we work in Euclidean signature but because the spacetimes we are interested in are static we can 
trivially Wick-rotate the time coordinate and obtain the solution in Lorentzian signature. 
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A = B = F = Qm$g we recover the Schwarzschild black string in D spacetime dimensions. 
To bring this latter solution to the standard form in Schwarzschild coordinates one performs 
the following change of coordinates: 




r, 



,£>-4 

o 



(3) 



r 



D-4 ' 



Therefore, the range of the y coordinate is y 6 [0, 1). The x coordinate will be our compact 
coordinate and we can take it to lie in the range 16 [0, 1] without loss of generality. In these 
coordinates y = is the horizon and y = 1 is the asymptotic region. Similarly, x = is 
the reflection plane and x = 1 corresponds to the periodic boundary. Therefore, with this 
particular choice of coordinates we have effectively set the length of the compact circle to be 



subject to suitable boundary conditions. Here £ a = g bc (T a bc — T a bc ) is the DeTurck vector, 
T is the Levi-Civita connection associated to our metric ([2]), and F is another connection 
associated to a reference metric g that we are free to prescribe]^] In our construction we have 
chosen g to be the metric of the Schwarzschild black string in D dimensions. 

With our choice of gauge, the boundary conditions that we impose are as follows. At the 
horizon (y = 0) smoothness of the metric requires that the functions X should be even in y 
and therefore we impose a Neumann boundary condition on all of them. At the reflection 
plane and the periodic boundary (x = 0, 1 respectively), all functions should be even. Finally, 
at infinity we impose that our solution is asymptotically KK which implies X = 0. With 
these boundary conditions the maximum principle of [T5J implies that £ a should vanish and 
all the solutions to the harmonic Einstein equations should in fact be Einstein (and more 
precisely Ricci flat in our case). Therefore, we can use £ a to monitor the numerical error of 
our solutions. 

To solve the equations numerically we have used a Chebyshev collocation approximation. 
For low dimensions and small A a single patch that covers the whole computational domain 
is enough to give accurate numerical results. However, as D increases it becomes harder to 
find accurate numerical solutions because the gradients of the functions X become larger and 
typically higher resolutions are required. To solve this problem we have used 2 conforming 
meshes or up to 3 non- conforming meshes to cover the domain]^] This has allowed us to 
get enough resolution in the regions where it is needed, which is basically where the string 
pinches in the A — > oo limit. For the results presented below the numerical error is estimated 
(monitoring £ a ) to be less than 1%. In addition, we have checked that our numerical solutions 
converge exponentially to the continuum which, given our numerical scheme, strongly suggests 
that there exist smooth NUBSs. 

The Einstein vacuum equations are scale invariant and we have exploited this freedom 
to fix the overall scale by setting L = 2. Then, in any given D, the family of NUBSs is 

5 In principle one could be more general and F need not be the Levi-Civita connection associated to g [7J. 
6 Ref. [6] followed a similar approach. 
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-- 2. 

With the ansatz pi) we solve (numerically) the harmonic Einstein equations [7] 



Rab — V( a ^fe) = 



(4) 
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parametrised by the inverse temperature (3 and therefore to move along the family we can 
vary this quantity^] With our ansatz for the metric (|2]), the inverse temperature of the solution 
is given by 

= 7 > (5) 

n — 1 

and therefore we can move along the branch of NUBS by varying r$. 

Finally we note that because in our construction the metric coefficients are even functions 
of y, we can analytically continue the range of the y coordinate to negative values. This 
corresponds to going across the bifurcation surface and to the other side of the Einstein- 
Rosen bridge with another asymptotically KK region being at y — > — 1. As we shall see in £j3] 
this will be useful in order to study the stability of the NUBSs via the local Penrose inequality. 



2.2 Results 

In this subsection we report on our results for NUBSs in D = 11, 12, 13, 14, 15 dimensions]^ 
As we have explained above, for fixed D, we keep the length of the asymptotic KK circle fixed 
and we move along the branch by varying the (inverse) temperature, starting at the GL point, 



2.2.1 Non-uniformity 

As just mentioned, we use (3 to parameterize our NUBS solutions when constructing them 
numerically. However, it is more usual to use A as the parameter and so we start in Figure [4] 
by showing the dependence of on A for D = 11, 12, 13, 14. These plots reveal the range of 
A for which we have constructed solutions. 

For D — 11, (3 is a monotonically increasing function of A, i.e., the temperature decreases 
with A. This is the same as happens in lower dimensions [3 El [7j. The plot flattens out 
beyond A ~ 1, suggesting that (3 is already close to the limiting value it attains at the merger 
point A = oo. 

For D = 12 we find qualitatively new behaviour. As we move along the branch of NUBSs 
starting from the GL point and increasing (3, we cannot find any solutions beyond (3 = (3q = 
0.6760 L. For this solution A = Ao = 0.258 and therefore it cannot be the end of the branch 
(which we expect to have A = oo). Instead it is a maximum of /3, i.e., a minimum of the 
temperature. Ref. [7j found a maximum /3 on the localised black branch for D = 5; our 
results show that this maximum of f3 switches to the NUBSs in D = 12j^] To find solutions 
with A > Ao we have used as initial guess for the Newton algorithm a NUBS with D — 13 
and a similar value of A EE 

7 Actually we'll have to be more careful for D = 12 because (3 does not uniquely characterise the solution. 

8 The results in D = 15 are qualitatively similar to those in D = 14 and we will not display them below. 
Our data extends up to A ~ 0.8. 

9 Here we assume that in the yet-to-be constructed localised black hole branch the area of the horizon 
behaves monotonically as a function of f3. 

10 For instance, a NUBS in D — 13 with A ~ 0.3 gives a good enough initial guess to find NUBSs in D = 12 
with A > Aq. 
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Figure 4: /3/L vs. A for D = 11,12,13,14. The overlapping curves correspond to data 
generated at different resolutions. In D — 11, (3 is a monotonically increasing function of A. 
This is the same behaviour as in the lower dimensional cases [5j|6l[7]. In D = 12, reaches 
a maximum at A = Ao = 0.258 and then it decreases, presumably up to the merger point. In 
D>13,/3/Lisa monotonically decreasing function of A. 



For D = 13,14,15 j3 is a monotonically decreasing function of A, i.e., the temperature 
increases with A. 
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2.2.2 Horizon area and mass 

In Figure [5] we plot Ah/L d ~ 2 vs. (3/L for various dimensions (Ah is the horizon area) 
and for completeness in Figure [6] we plot M/L D ~ 3 vs. (3/L. The first law implies that the Ah 
and M exhibit the same qualitative behaviour. Using these data, we have checked that the 
error in the Smarr relation is consistent with our estimates based on the norm of the DeTurck 
vector £ a . 

For _D = 11, Ah is a monotonically increasing function of (5 (or A), at least up to A ~ 
2.4 but it is reasonable to expect that this is true for all values A up to the merger point. 
This qualitative behaviour is the same as in the lower dimensional cases discussed previously 

For D = 12 we find qualitatively new behaviour. As we move along the branch of NUBSs 
starting from the GL point and increasing (3, the area of the horizon increases up to, and 
beyond, the minimum temperature (maximum 0) solution. However, it reaches a maximum 
at (3 = f3\ = 0.6757 L. The value of A at this point is A = Ai = 0.332. The horizon area 
decreases for larger A. By the first law, the mass also has a maximum at this value of (or A). 
As in the localised black hole case for D = 5 [7], (3\ is close to (3 , but they do not coincide: 
ft/A, ~ 0.9995 ± 0.0001. 

For D = 13 we find that there is again a maximum area/mass solution, which has /3 = 
Px = 0.6356 L and A = Ai = 0.135. 

For D = 14,15, the area/mass decreases monotonically with A, presumably all the way 
to the merger point. This was the expected behaviour and our results confirm the expec- 
tations well in the non-perturbative regime. However, as D increases it becomes harder to 
find accurate numerical solutions (within reasonable time and with reasonable computational 
resources) and we have only been able to explore the branch of NUBSs up to A ~ 1. 

2.2.3 Comparison with uniform string 

In Figure [7] we plot the difference in horizon area between NUBSs and uniform black 
strings with the same mass for various space-time dimensions. For D — 11 we find that the 
NUBSs always have lower horizon area than uniform black strings with the same mass, as in 
the previously studied lower dimensional cases |6j [7j. This implies that a D < 11 NUBS 
cannot be the endpoint of the GL instability, as explicitly demonstrated for D = 5 in Ref. 

For D = 12, 13 we find that there exists A2 > Ai such that for any A > A2, the NUBSs have 
greater entropy than the uniform black strings with the same mass. We find A2 = 0.733 in 
D = 12 and A2 = 0.192 in D = 13. Therefore, it is natural to expect that NUBSs with A > A2 
could be the endpoint of the GL instability, at least in some region of the phase diagram. 

For D = 14, 15, we find that any NUBS has greater entropy than a uniform black string 
with the same mass, which agrees with the perturbative results of Ref. [I] and extends them 
well into the non-perturbative regime. 
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Figure 5: Ah/L d ~ 2 vs. j3/L for D = 11, 12, 13, 14. The black curves are the uniform strings, 
the red curves are our results for NUBSs. For D = 11 the area of the horizon is an increasing 
function of (3, just as in lower dimensional cases previously considered [3, El [7]. For D = 12 
there exists a maximum inverse temperature at /3q/L ~ 0.6760 and a maximum of the area 
at fli/L ~ 0.6757. In D = 13 the maximum of the area is at Pi/L ~ 0.6356 and in D > 14 
the area always decreases as (3 decreases. By the First Law the mass has the same qualitative 
behaviour as the horizon 
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Figure 6: M/L D ~ 3 vs. (3/L for D = 11,12,13,14. The black curves are the uniform black 
strings. The behaviour of M as a function of (3 in various dimensions is similar to that of 
horizon area Ah, as dictated by the first law. 

2.2.4 Heat capacity 

Let cl = dMj dT be the heat capacity, with the subscript L referring to the fact that we 
are keeping the asymptotic length of the KK circle fixed. 

For D = 11, the horizon area increases with (3 which implies that cl is negative. This is 
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Figure 7: Comparing the horizon area of a NUBS with that of a uniform string of the same 
mass. The vertical axis is Ah/A° h — 1 where Ah is the NUBS horizon area and A° H the 
horizon area of a uniform string with the same mass. The overlapping curves in this figure 
correspond to data at different resolutions. In D = 11 we observe the same behaviour as in 
lower dimensions: the NUBSs have lower entropy than the uniform black strings with the 
same mass. In D = 12, 13 there exists a A = A2 such that for A > A2 the NUBSs have greater 
entropy than uniform black strings with the same mass. For D > 14 and any A > the 
NUBSs always dominate the microcanonical ensemble. 
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the same as found previously for D < 11 dimensions [SJ El E] • 

For D = 12, writing cl = (dM / dX) / (dT / dX) shows that cl changes sign by diverging at 
the minimum temperature solution A = Ao and changes sign by vanishing at the maximum 
mass solution A = Ai. Hence cl is negative for A < Ao, diverges at A = Ao, is positive in the 
small range A < A < A 1; vanishes at A = Ai, and is negative for A > Ai, presumably all the 
way to the merger point and also along the localised black hole branch. Ref. [7] found that for 
D = 5 there is a portion of the localised black hole branch which has positive heat capacity. 
Our results suggest that solutions with positive heat capacity remains on the localised black 
hole side up to D = 11. However, for D = 12 they switch to the NUBS side. 

For D = 13, cl is positive for A < Ai, vanishes at A = Ai, and is negative for A > Ai. 
Presumably the localised BHs have negative c^. 

For D = 14, 15, all NUBSs that we have been able to construct have negative heat capacity. 
Presumably this is true also of the localised BHs. 

2.2.5 Free energy 

In Figure [8] we plot, for various D, AF = F — F° where F is the free energy of a NUBS 
and F° the free energy of a uniform black string with the same temperature. 

For D — 11, the NUBSs always have greater free energy and therefore they never dominate 
the canonical ensemble. 

For D = 12 the free energy exhibits a cuspy behaviour. Small A NUBSs have greater free 
energy than uniform black strings. As A increases, AF reaches a maximum at the minimum 
temperature A = Ao, i.e., /3 — /3q — 0.6760 For A > Ao, AF decreases, eventually 
becoming negative. It becomes zero at /3 = 0.6748 L, so NUBSs with smaller /3 are preferred 
over uniform strings in the canonical ensemble. 

For D = 13, 14, 15 the NUBSs always have lower free energy than the uniform strings. 

2.2.6 Negative modes 

Next we consider the spectrum of the Lichnerowicz operator A^ restricted to modes which 
respect the U(1) T x SO(D — 2) isometries of the background. We will only be interested in 
the negative modes of A^, which are an invariant of the manifold in question. Our results for 
different D are displayed in Figure [9j 

For D — 11 we find that NUBSs have two negative modes, at least for all values of A that 
we have been able to explore. The mode with the most negative eigenvalue is continuously 
connected (as A — > 0) to the uniform black string negative mode arising from the Euclidean 
Schwarzschild negative mode. The other NUBS negative mode reduces to the GL zero mode 
as A 0. It is interesting to note that there are no other zero modes along the NUBS 
branch, just as in D = 5 [TJ. However, we do not see that the most negative mode becomes 
increasingly negative along the branch; in fact, both negative modes appear to remain finite 

11 To see why these extrema must coincide, note that dF/dX = dM/dX — TdS/dX — SdT/dX. The first law 
implies that this vanishes at A = Ao since dT/dX = 0. Hence F is extremized at A = Ao- Furthermore, F° is 
a single-valued function of T so dF a /dX = (dF° / dT)dT / dX, which also vanishes at A = Ao- Hence, viewed as 
a function of A, F° also is extremized at Aq. Therefore AF also is extremized. 
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Figure 8: Difference between the free energy of NUBSs and uniform black strings for the same 
temperature. The multiple curves correspond to data at different resolutions. In D = 11 
NUBSs always have greater free energy than uniform strings with the same temperature. 
In D = 12, AF reaches a maximum at the minimum temperature and it becomes negative 
for (3 < 0.6748 L, indicating that for these temperatures the NUBSs dominate the canonical 
ensemble. In D > 13 the NUBS always have a lower free energy than the uniform black 
strings with the same temperature. 
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Figure 9: Negative modes of the Lichnerowicz operator in various dimensions. For D — 11 
NUBSs have two negative modes all along the branch. For D = 12, weakly non-uniform strings 
(i.e., small A) have two negative modes; at the minimum temperature (ft = /3 max ) the mode 
with the least negative eigenvalue turns into a zero mode. For larger values of A there is only 
one negative mode. For D > 13 NUBSs have only one negative mode. 



as A increases. Of course, we have only been able to construct solutions up to moderate values 
of A and it is entirely plausible that the mode with the most negative eigenvalue diverges at 
the merger point, as in D = 5 [7]. 
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Figure 10: Embedding of the horizon geometry of a D = 11 NUBS with A = 2.41 into M . 
On the left figure we have plotted the radius vs. X as a surface of revolution. On the right 
we compare the embedded geometry with that of a Ricci flat cone. The geometry starts 
approximating the Ricci flat cone metric of Ref. [8]. 



For D = 12 we find that NUBSs with A < Ao have two negative modes, exactly as in 
the lower dimensional cases. However, one of the negative modes reduces to a zero mode at 
A = Ao and is absent for A > Ao- The zero mode at A = Ao is obtained by an infinitesimal 
variation of parameters of the NUBS solution. Only at an extremum of temperature does such 
a variation of parameters respect the boundary conditions imposed on a Euclidean negative 
mode, namely that it should preserve the temperature. Hence a Euclidean negative mode can 
reduce to a zero mode corresponding to a variation of parameters in the background solution 
only at an extremum of temperature [T7] . For A > Ao we find that there is only one negative 
mode and no other zero modes. 

For D > 13, we find that all NUBSs have one, and only one, negative mode. This mode is 
continuously connected (as A — > 0) to the negative mode of the uniform black string inherited 
from the Schwarzschild solution. 



2.2.7 Local geometry 

Finally we consider the local geometry near the minimum radius S D ~ 3 of the horizon. 
Recall that Ref. [8] predicted that precisely at the merger point, the local geometry should be 
given by that of a cone over S 2 x S D ~ 3 , where first sphere comes fibering the U(1) T and the 
second sphere comes from the axisymmetry of the problem. Therefore, one would expect that 
the local geometry near the minimal sphere of a NUBSs approximates that of a cone only if 
the non-uniformity parameter A is sufficiently large. Following Refs. [18j |6], in Figure 10 we 



16 



have embedded the geometry of the horizon of a D = 11 NUBS with A = 2.41 into Mr 2 

dsln-2 = dX 2 + dY 2 + Y 2 dQ 2 D _ 3 , (6) 

and projecting onto an equatorial circle of the S D ~ 3 in the plots. It is interesting to note 
that as D increases the geometry of a horizon of given A looks "rounder" (compare with 
the embedding plots in Refs. [5j [6]). To compare with the cone geometry we have also 
embedded the latter; as Figure [TO] shows, the geometry of the horizon of the A = 2.41 NUBSs 
is reasonable close to that of the cone, but given that A is not large one should not expect 
very good agreement. For the solutions that we have constructed with D > 11, A is even 
smaller and the geometry near the minimal horizon sphere is still quite far from the cone. 



3 Local Penrose inequality 

3.1 Introduction 

To investigate stability, it is convenient to label NUBSs by their mass M instead of A. 
For those D in which the mass exhibits a maximum at A = Ai, we divide NUBSs into two 
families: "thin" non-uniform strings with < A < Ai, and "fat" non- uniform strings with 
A > Ai. Solutions belonging to each family are uniquely labelled by M. Stability requires 
that an arbitrarily small perturbation leads to an arbitrarily small change in the solution, so 
a stable thin (fat) solution cannot evolve to a fat (thin) solution when perturbed. Hence, in 
investigating stability we can regard thin and fat strings as distinct families of solutions. Let 
y 4-NUBs(-^0 denote the horizon area of the particular family of interest. 

We will study stability using the method of local Penrose inequalities introduced in Ref. 
|13j . Consider a 1-parameter family of time-symmetric initial data for the Einstein equation, 
specified by a parameter e, with the following properties: (a) it is asymptotically flat in the 
KK sense, with the KK circle having circumference L at infinity; (b) for e = it reduces to the 
initial data on a constant r slice (i.e. the Einstein-Rosen bridge) of a NUBS with mass Mj. 
For small e, our initial data describes a perturbation of the NUBS with mass Mj. Assume that 
this solution is stable. This means that, for sufficiently small e, the perturbation will disperse 
through radiation to infinity and across the horizon, and the spacetime will settle down to 
a new stationary solution, with mass M/(e), close to Mj. In general, the final state might 
have some angular or linear momentum and would belong to a stationary generalization of 
the static NUBS family of solutions. Since we know nothing about this generalization, we can 
exclude this possibility by imposing the addition restrictions: (c) the initial data preserves the 
SO(D — 2) symmetry of our NUBS solutions; (d) the initial data has a reflection symmetry 
around the KK circle. Our final state must then belong to the same family (thin or fat) of 
static NUBS solutions as our initial solution. 

Let Aapp(e) denote the area of the apparent horizon of the initial data, and Ai(e) the area 
of the intersection of the event horizon with the initial data surface. The apparent horizon lies 



within the event horizon, The time symmetry of our initial data implies that the apparent 



12 This assumes that cosmic censorship is not violated in the evolution, which is justified because a violation 
of cosmic censorship would be an instability, contradicting our assumption of stability. 
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horizon is a stable minimal surface and hence has smaller area than any surface which encloses 
it. Hence 

A app (e) < Ai{e) < A NVBS (M f (e)) < A NUBS (M( e )) (7) 

The second inequality follows from the Hawking area theorem. In the third inequality, we use 
Mf(e) < M(e) where M(e) is the mass of our initial data. This follows from the decrease of 
Bondi energy (i.e. gravitational waves carry away energy). The above inequality is referred 
to as the local Penrose inequality for the NUBS. The word "local" refers to the fact that it 
should hold only for arbitrarily small e. It is a necessary condition for stability of the NUBS. 

This inequality is saturated by initial data corresponding to the NUBS solution. Hence a 
stable NUBS is a local maximum of horizon area in the space of initial data of given mass, 



Q > 0, Q = M(0) - 7 Ti app (0) - — M(0) 2 



and satisfying the above conditions, 13 The first law implies that an unstable NUBS is an 
extremum of horizon area at fixed mass, but this extremum is not a local maximum. 

We are interested in arbitrarily small e so we expand both sides of this inequality in e. At 
zeroth order we have M(0) = Mi, y4 app (0) = Anubs(M0 an d so th e inequality is saturated. 
At first order the inequality remains saturated because of the first law of black hole mechanics 
M(0) = (T/4)v4 app (0) where a dot denotes a derivative with respect to e and T is the Hawking 
temperature of the NUBS with mass Mj. At second order in e, the local Penrose inequality 
reduces to 

1 

where cl = (dM/dT)L denotes the heat capacity of the NUBS (with L held fixed). 

To summarize, a necessary condition for stability of the NUBS is that initial data de- 
scribing a small perturbation of the NUBS must satisfy the inequality Q. The definition 
of Q involves quantities which depend on the perturbation to second order in e. However, 
the constraint equations governing second order perturbations of initial data reveal that such 
perturbations are sourced by quantities quadratic in first order perturbations and it can be 
argued using the first law that Q itself depends only on the first order perturbation |13j . 

Is the condition ^ also sufficient for stability? Ref. [19] argued that it is both necessary 
and sufficient for the stability of linearized perturbations. This was done by defining a quantity 
£ that can be expressed as an integral over a spacelike surface of a quantity quadratic in 
a linearized perturbation, and showing that £ is a non-increasing function of time. If all 
perturbations have £ > then linearized perturbations cannot grow, strongly suggesting 
stability. However, if there exists a perturbation with £ < then such a perturbation cannot 
decay, strongly suggesting instability. Ref. [19] proved that there exists a perturbation for 
which £ < if , and only if, there exists a perturbation for which Q < 0. 

Let's apply the above arguments to a NUBS solution with D < 13. Consider first the 
solution corresponding to the intersection of the uniform and non-uniform families of solutions. 
We can regard this as a NUBS with A = 0. Consider the perturbation of this solution 
corresponding to an increase in mass, moving along the uniform string branch. This gives 

13 Conditions (c) and (d) are not required to reach this conclusion: they can be eliminated if one knows 
that a generalization of the static NUBS with non-vanishing linear or angular momentum has lower horizon 
area than a static NUBS of the same mass. This follows from the first law and the assumption that the linear 
and angular velocities of the horizon have the same sign as the linear and angular momentum respectively. 
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a uniform string with a greater horizon area than a NUBS with the same mass. Hence this 
perturbation violates the local Penrose inequality for the NUBS family. So if we regard this 
solution as a member of the NUBS family then it is unstable. Now consider a NUBS with very 
small non-zero A. Construct initial data describing a perturbation of such a solution, such 
that it reduces to the perturbation just discussed in the limit A — > 0. Expanding Q in powers 



of A givet Q = Q + XQi + A Q2 + We've just seen that Q < 0. Hence for sufficiently 



small positive A we will have Q < 0. Therefore a D < 13 NUBS with small positive A will be 
unstable. 

We have seen that, for D = 12, 13, there exists a maximum mass along the branch of 
NUBS solutions. The same occurs in other contexts, e.g. rotating neutron stars or black 
rings at fixed angular momentum J. In these other systems, a heuristic "turning point" 
argument suggests that an unstable mode appears when one passes through the maximum of 
mass, on the side in which the heat capacity is positive [E1EQ]. The local Penrose inequality 
provides a more rigorous justification for this argument |13j . At a maximum of mass cl (cj 
in the case of rotating solutions) changes sign by passing through zero. Hence the final term 
of ^ diverges, and is negative on the side in which cl is positive. This strongly suggests 
that solutions close to the maximum mass solution, with positive cl should be unstable. Note 
that it does not suggest any change in stability at the minimum temperature solutions, where 
cl changes sign by passing through infinity (this is why there is no instability of the Kerr 
solution when the heat capacity at fixed angular momentum changes sign [21]). 

Our D = 12, 13 NUBS solutions exhibit a maximum mass at A = Ai, with cl positive for 
A < Ai, so the above argument strongly suggests instability for A close to, but less than, Ai. 



3.2 Instability of non-uniform black strings 

We have argued above that NUBSs are unstable if D < 13 and A is infinitesimal, or 
if D = 12, 13 and A is slightly less than Ai. In this section we will use the local Penrose 
inequality to investigate the stability of NUBSs with other values of D and A. 



3.2.1 Simple initial data 

Let h ab the induced metric on a surface of constant r in a NUBS solution. The extrinsic 
curvature of such a surface vanishes, so we consider time-symmetric initial data and hence 
the momentum constraint is trivially satisfied. Construct new time symmetric initial data, 
describing a perturbation of this solution, by conformally rescaling the induced metric: 

h ab = ^ D ~^h ab , K ab = 0. (9) 

The Hamiltonian constraint then becomes 

V 2 ^ = (10) 

where V is the covariant derivative with respect to h ab . Expanding \1/ in terms of a small 
parameter e as \I/ = 1 + e ^(0) + e 2 ^(0)/2 . . .J^Jwe can construct, perturbatively, initial data 

14 One can show Q\ — but that doesn't matter here. 

15 Here the dots ' denote derivatives with respect to the parameter e. 
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by solving 

V 2 ^ = V 2 ^ = . (11) 

As explained in Ref. p2], the homogeneous solution for \1/ does not contribute to the 
local Penrose inequality. Therefore we choose the trivial solution ^ = 0. For the first order 
solution we demand that the resulting induced metric is asymptotically KK, which amounts 
to impose \I> = (9(l/r n_1 ) = 0(l — y) near infinity. With this decay at infinity, the asymptotic 
length of the KK in the perturbed space-time is the same as that of the NUBS background. 
In addition, we place the inner boundary at y = yt < and impose Dirichlet boundary 
conditions there, i.e., we specify ^\ y=yb . We can do this because in our numerical solutions 
for the background NUBSs all functions are even in y so we can analytically continue the 
range of y to negative values. Then, y = is the bifurcation surface and the region of the 
space-time with y < corresponds to another asymptotically KK region. Finally, we impose 
reflecting boundary conditions on \I/ on both the reflection plane (x = 0) and on the periodic 
boundary (x = 1). Solving the equation for \I/ with these boundary conditions, we can obtain 
the variation of mass, the variation of the area of apparent horizon and finally Q. We describe 
the details of the calculation of Q in Appendix [A} 

For the results presented in Fig. ITT] we have imposed the simple boundary condition 



ty\ y=yb = 1, and this turns out to be sufficient to detect instabilities for all values of A 
that we could explore. We also considered more general boundary conditions ^\ y=yb = 1 + 
Y2n=i Cn cos ( n7rx ) i n or der to look for instabilities for those values of A for which Q > 
with the previous simple boundary condition. Here the integer N represents a truncation of 
the Fourier expansion and we fixed the normalization of \I> by choosing the first term in the 
expansion to be unity. To look for initial data which satisfies Q < 0, we seek to determine 
the {c n } which minimize the dimensionless quantity 

Q = -^, D 2 = £ /2 dx*(y = 0,xy . (12) 
We varied the number of Fourier coefficients N in the previous expansions, but the qualitative 



results presented in Fig. 11 do not change. Finally, we put the inner boundary at y& = 
— 1.25 x 1CT 3 , but have checked that our results do not depend on the specific choice of y^. 

To proceed with the calculation of Q, we first interpolate the (numerical) background solu- 
tion using a Lagrange polynomial and on top of this background we consider the perturbations 
as we described above. In this calculation, it is crucial to ensure that the numerical error of 
the background does not affect the value of Q. To do so, recall from the discussion in ^2] that 
<f) = £ a £ a provides a measure of the numerical error of the background. Therefore, we have to 
estimate the value of <fi = ^acceptable such that for (f> < (^acceptable the value of Q is independent 
of cf). We have made this estimate by comparing the values of Q obtained for backgrounds 
computed at two different resolutions (and hence different values of 4>) and we found that for 
0max < 10~ 3 the value of Q changes by less than 10% as the resolution of the background 
changes. This estimate is independent of the number of space-time dimensions. Here </> max 
denotes the maximum value of <fi over the whole computational domain. The calculation of 
Q has its own numerical error, which we estimate by monitoring 5 = \M — 7M app /4|/(.DM), 
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which should vanish by the first law of black hole mechanics. For all results presented the 
error in the background turns out to be the dominant one and 6 < 0.01 in all cases. 



Our results are presented in Figure 11 where we plot Q against the non-uniformity pa- 
rameter A in various dimensions. For D = 11, we see that Q is negative for any value of A 
for which we have acceptable numerical errors. Thus all D = 11 NUBSs we have checked 
are unstable, even well in the non-perturbative regime (i.e., large A). This is one of the main 



results of this paper. _ We have also checked that NUBS in D = 5, 6, 8, 10 are also unstable 
for different values of A, which are not necessarily small. Therefore, it is reasonable to expect 
that all non- uniform black strings in 5 < D < 11 are dynamically unstable. 

For D = 12, 13, we find that Q is negative for A < Ai, which implies that all "thin" 
non-uniform strings in D = 12, 13 are unstable. Note the divergence in Q at A = Ai, where 
cl vanishes p 7 ] Thus our results confirm the existence of an instability when A is close to, 
but less than, Ai. More importantly, they demonstrate that an instability is present for all 
A < Ai. For A > Ai, we minimised Q as explained above (with N = 3,4) and we always 
found that Q is small and positive and hence, there is no indication of an instability in "fat" 
non-uniform strings. Given the ease with which we found an instability for thin strings, this 
strongly suggests that fat black strings are stable. 

Finally, for D = 14 we found that after minimisation Q is positive for any A, which again 
strongly suggests that these solutions are stable. We found the same result for D = 15. 

3.2.2 More general initial data 

We have also considered a more general class of time-symmetric initial data defined as 
follows. Instead of taking h a b to be the induced metric on a surface of constant r in the 
NUBS spacetime, consider the deformed metric 



2 S 

hob -> h ab dx a dx b = rpe 2 e" 2,I>/(D - 3) dQ 2 D _ 3) + e* 



e A dx 2 + 4r °^ e 3) (dy + yf(y) F dxf 

f{y) D ~ 4 



(13) 

where the function <§>{x,y) describes the deformation: for $ = the metric (13) reduces to 
the induced metric of a NUBS. Note that for an arbitrary $ the induced Ricci scalar R is 
non-vanishing because the Hamiltonian constraint is not satisfied. 

To construct initial data that satisfies the constraints we proceed as before by conformally 
rescaling the induced metric as in rt9|) and solving the Hamiltonian constraint in terms of the 
conformal factor: 

V 2 fr= ?~ 3 ^ RV, (14) 
4(D-2) K J 

where V is the covariant derivative with respect to h a b- Note that $ is a free function; once 



we have chosen $ and specified boundary conditions, eq.(14) uniquely determines \&. The 



idea now is to expand both $ and $ in terms of the small parameter e. So as before we have 



16 Ref. [TT] argued that slightly NUBS should be unstable by continuity. 

17 A similar situation was found in Ref. [T5] for black rings: "fat" rings were found to be unstable and "thin" 
rings were found to be stable and Q goes through infinity when moving from one branch to the other because 
the specific heat at constant angular momentum vanishes at the maximum area (at fixed mass) solution. 
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(a) D = 11 



(b) D = 12 




(c) D = 13 



(d) D = 14 



Figure 11: Q vs. A for D = 11, 12, 13, 14. In Figures (b) and (c) (D = 12, 13 respectively) 
we have included insets showing a zoom in of the region with A > Ai ("fat" NUBSs). For 
D = 11, Q is always negative, which implies these that NUBSs are unstable. For D = 12 and 
13, Q is negative for A < Ai and positive for A > Ai. Thus, "thin" NUBSs are unstable, but 
for "fat" NUBSs, there is no indication of instability. For D = 14, Q is positive for any A and 
there is no indication of instability. 



if/ = 1 + e *(0) + e 2 *(0) /2 + . . ., and we demand that $ = 0(e) so that in the e -)> limit we 
recover the NUBS. For simplicity (and because we can freely specify this function) we choose 



$ = e$, and all the higher order terms in e identically vanishing. Then, expanding (14) in 
powers of e up to second order gives 

V 2 ^ = D ~ 3 A 

4(D r 2) . as) 

v 2 ^ = ^ _ 3 2) [& + 2 k (* + $)] + 2 /^(<9» (a 6 *) . 
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Solving (15) subject to boundary conditions as in the previous subsubsection we have 
constr ucted, perturbatively, initial data for different choices of $. We have proceeded as in 
£ 3.2.1 and expanded the inner boundary condition for \1/ in a Fourier series and then minimised 

y=Vb = (which we have), 
Therefore, in general we have M(0) ^ 0. 



Q. Note that now \1/ is sourced by $ and so even if we impose \l/ 
this function will be non-vanishing 



The results that we have obtained for the Penrose inequality in various dimensions us- 



ing this more general class of initial do not qualitively differ from those reported in £3.2.1 



Therefore we do not display them here. It is worth emphasising that no matter the choice of 
$, we always find that Q is positive in D = 12, 13 for "fat" NUBSs and in D = 14, 15 for all 
NUBSs. This is further evidence that such solutions ought to be stable. 
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Details of the calculation of the local Penrose in- 
equality 

In this appendix we give more details on the calculation of Q (see Eq. rt8| for the definition) 



for the class of initial data in £ 3.2.1 To do so, we need to calculate the variation of the mass 
and area of apparent horizon induced by the perturbation. 

Lets first recall the expressions for the mass in asymptotic Kaluza-Klein spacetimes, i.e., 
the metric at infinity approaches Mo-i x S 1 , where Mu-x denotes the (D — l)-dimensional 
Minkowski space. In the linearised regime one finds that the spatial components of the metric, 
in a Cartesian coordinate system, are given by [22]: 



where r = y ^2f =1 2 (x 1 ) 2 is the Euclidean radial distance, z ~ z + L is the compact direction 
and c t and c z are two constants. Then the mass is given by 

M= I ^[(D-3)c t -c z ], (17) 

18 As explained in Ref. |13] , changing the inner boundary condition for ^ amounts to modify the homoge- 
neous solution to the equation which, by the first law, cannot affect the Penrose inequality. 
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where £Ids is the area of a (D — 3)-dimensional sphere. Given the form of our initial data, 
Eq.([9|, we can easily calculate the variation of these parameters due to the perturbation 

c z = ^^(r D ~ 4 n. =00 , c t = (r D -^) r=00 . (18) 

Since we set ^ to be zero, we have c z = <it = 0. Thus, there is no second order perturbation 
in variation of the mass. The first order variation of mass due to our perturbations is 

M = (D ~ 3)2 ~ V - 4 ^ . (19) 

4tt D-3 V ; V ; 

Now, we consider the variation of area of apparent horizon. In Ref.|13j. it was shown 

that the apparent horizon coincides with the minimal surface for axisymmetric perturbations. 

Thus, we calculate area of the minimal surface instead of that of the apparent horizon. We 

denote the location of the minimal surface as y — Y(x). Then, the induced metric on the 

minimal surface is 



-y IJ dx I dx J = ^ D -V 



-^r- dnl_, + e A dx 2 + 4r ° 2 tf 3) (Y> + Yf{Y)F) 2 dx 2 



(20) 



where ' = d/dx. In the above expression, the functions A, B, F, S and \I/ are evaluated at 
y = Y(x) and they depend only on x. The area of the minimal surface is given by A app = 

J d D ~ 2 Xy/j where 7 = det(7/j). We expand Y(x) and v4 app in term of e as Y(x) = eY(x) H 

and y4 app = Abg + e^4 apP + e 2 ^4 apP /2 + • • • where Abg is the area of the event horizon in the 
background geometry. Then, we obtain first and second order perturbation of the area of the 
minimal surface as 

i apP = ^f^-a rt 3 J dx e^ D ~ 3 ^ . (21) 
i apP = JW?" 3 [ dx e ^ D ~ 3 ^ 2 Ur 2 Ae B - A (Y' + YFf + 2( ^~ 3) F 2 

Jh L d — a 

1 (22) 



where functions A, B, F, S, ^> and their derivatives are evaluated at background horizon y = 0. 
Since Y(x) extremizes the area functional A app , we have the equation for Y(x) as 

- g r 2A e -(A+(D-3)S)/2^ e (-A+2B+(D~3)S)/2^ Y , + y F)']' 



+ 8r 2 Ae B - A {Y' + YF)F + {d 2 A + (D — 3)d 2 S)Y + K ' djb = . (23) 

1 J O 

Using this equation, we obtain a simple expression for v4 app as 

K w = ^ffWa J dx e ^ D -^' 2 [Yd y i> + §^ 2 ] • (24) 

In the calculation of Q we do not need the first order variation of the area of the ap- 
parent horizon A app . However, we use this quantity to check the first law of the black hole 
thermodynamics M — TA &PP /A = 0, which provides an estimate of the numerical error in our 
calculations. 
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